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Abstract. Inelastic and diffractive production of J/-0 mesons at HERA is reviewed. 
The data on inelastic photoproduction are described well within errors by the Colour 
Singlet Model in next-to-leading order. A search for colour octet processes predicted 
within the NRQCD/factorisation approach is conducted in many regions of phase 
space. No unambiguous evidence has been found to date. Diffractive elastic production 
of J/V' mesons has been measured in the limit of photoproduction (Q^ ~ 0) up to the 
highest photon proton center of mass energies. The increase of the cross section is 
described by pQCD models. At larger Q^, 2 < < 100 GeV^, the W^p dependence 
is found to be similar to that observed in photoproduction. First analyses of data at 
high t, \t\ < 20 GeV^, yield a powerlike dependence on \t\. A LO BFKL calculation 
gives a good description of the data. 
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1. Introduction 

The main interest in leptoproduction of J/ip mesons is a clarification of the production 
process. In the inelastic regime, many models have been proposed, and predictions 
within the theoretically favoured NRQCD/factorisation approach could not yet be 
identified nor could they be ruled out. In the diffractive regime, models based on pQCD 
have proven to work in principle and more detailed studies are underway experimentally 
and theoretically. 

This review will cover J /if) production at HERA in inelastic and diffractive 
processes. Inelastic J/ip production at HERA has been studied in the limit of 
photoproduction (Q^ — > 0) and for > 2 GeV^. Data at high Q"^ were published by 
the HI collaboration [|l| and new theoretical calculations are discussed in L. Zwirner's 
contribution [Q. Here I concentrate on photoproduction. 

Diffractive elastic J/ip production is studied with large statistics in the limit of 
photoproduction and also at > 2 GeV^, where the W^p and dependences are 
analysed. Proton dissociative data at high values of t, the momentum transfer squared 
at the proton vertex, are also analysed. 

2. Inelastic Production of J /if) Mesons 

The data from the HI and ZEUS collaborations on inelastic J / ip photoproduction 
are based on integrated luminosities of 22 pb~^ and 38 pb~\ respectively, collected in the 
years 1996/97. The kinematic regions are characterized by < 1 GeV^, 50 < W^p < 
180 GeV, 0.4 < 2 < 0.9 for the ZEUS data and 60 < W^p < 180 GeV^ 0.3 < ^ < 0.9 
for HI. The usual kinematic quantities = —{k — k'Y and W^^ = + P)"^ used, 
where k, k', q and P denote the four-momenta of the incoming and outgoing electron, 
of the exchanged photon and the incoming proton, respectively. An important variable 
is 

q-P^^ 

z = 




Figure 1. Generic diagrams for J/ip production mechanisms at HERA: a) diffractive 
elastic process with two gluon exchange. Inelastic processes via boson gluon fusion 
with b) direct and c) resolved photons. 
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where is the four- momentum of the J /if) meson. In the proton rest frame z is the 
relative energy of the J /ip with respect to the photon, z k, The variable z is used to 

7 

distinguish diffractive processes, where z ~ 1, from inelastic processes, where z < 0.9. 

The inelastic production process at HERA is dominated by boson gluon fusion, 
see Figure |l]b,c. At medium z, direct processes dominate, where the photon couples 
directly to the charm quark, while at small z resolved photon processes are expected to 
contribute, where the photon interacts via its hadronic component. 

The theoretical description of J / ip meson production has undergone a change in the 
last few years. The Colour Singlet Model (CSM) has been replaced by an effective 
field theory approach: non-relativistic QCD (NRQCD) and factorisation ||^. This 
approach contains the CSM but predicts additional contributions. In the CSM, the 
cc pair is produced in the hard scattering with the observed quantum numbers of the 
J/ip meson. The colour singlet state is achieved by emission of an additional gluon, i.e. 
the process is •yg ^ g + cc [1,^ In the NRQCD/factorisation approach the cc pair 
can also be in a colour octet state where different angular momentum states are also 
allowed. 

At HERA we are in a unique position since for photoproduction full next-to- 
leading order calculations are available in the colour singlet model, which reduces the 

X The '1' or '8' denote the colour state of the cc pair and spectroscopic notation is used for the angular 
momentum, '^^'^^Lj. 
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Figure 2. Cross sections as functions of W-yp, and pt,^^. The data from HI |^ and 
ZEUS H, Q are shown with NLO CSM calculations. The band reflects the uncertainty 
due to the charm mass (1.3 < < 1.5 GeV) and 0.1175 < as{Mz) < 0.1225. The 
renormalisation and factorisation scales are fi^ — 2m^. (Figures from m.) 
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uncertainties in the normalisation of the predictions. So we will first show a comparison 
of the data with these calculations before we address possible contributions of colour 
octet (CO) intermediate states. 

2.1. Comparison of data to Colour Singlet calculations 

In Figure ^,b the z and the W^p distribution are shown for data from ZEUS 0] and 
HI which agree well with each other. The prediction of the CSM in next-to-leading 
order ||^ is overlaid. The band reflects the major theoretical uncertainties as detailed 
in the figure caption. A good description of the data is achieved. The necessity of 
including NLO contributions is demonstrated in the distribution of pj^^ in Figure |^c. 
The theoretical band including NLO contributions describes the data well and is, at 
^ Pt,t(> ~ 4 GeV, roughly a factor of 10 above the LO prediction. The discrepancy 
found in pp ^ J/ip + X at values of ~ 20 — 30 GeV between the data and the 
LO CSM calculation which led to the introduction of CO contributions within the 
NRQCD/factorisation approach was approximately a factor ~ 40 0. 

A different approach to describe the data within the CSM is shown in Figure ^. 
Here the HI data are shown with the results of two Monte Carlo implementations 
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Figure 3. Differential cross sections for photoproduction of J/tp mesons as functions 
of z, Wjp and and the ratio of Data/MC for the p^ bins, normahzed in the third 
bin. The predictions are from the CCFM Monte Carlo generator CASCADE and the 
Colour Singlet Model EPJPSI (the latter scaled by a factor 1.7) are compared to HI 
data with pt > 1 GeVand z < 0.8. 
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of the Colour Singlet Model in leading order. The dashed histograms represent the 
prediction from the EPJPSI program which uses "standard" parton distributions 
integrated over the transverse momentum and evolved with DGLAP equations. The 



program CASCADE |]rT| uses the CCFM evolution equation and "unintegrated" parton 
distributions. These have a distribution in transverse momentum, fcf, and are obtained 
from a fit to HI data on F^. The J /ip data in Figure |^(note that now a cut z < 0.8 has 
been applied) are in general well described by the CASCADE model except at very large 
z values, where relativistic corrections are missing. The pf^ distribution of CASCADE 
agrees better with the data than EPJPSI, but there is still a systematic discrepancy 
towards high pt. 



2.2. Search for Colour Octet Contributions 

In Figure ^ the same differential cross sections, da/dz from the HI and ZEUS data 
as before, are shown again as function of z. The theoretical predictions within the 
NRQCD/factorisation approach in leading order are shown in comparison. These 
include the direct photon contributions as well as resolved processes. At z > 0.3, 
direct photon contributions dominate and the CO contribution is dominant at 2; > 
0.6. The uncertainties of the summed prediction, shown as a band, are mainly due 
to the uncertainties of the long range matrix elements, the parameters which, in the 
NRQCD/factorisation approach describe the transition of the cc state to a physical J/ifj 
meson. 
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Figure 4. Differential cross section for photoproduction of J/tp mesons as a function 
of z. In a) the predictions from Kraemer Q are shown, in b) those from Kniehl et 
al. [ p^ . Both predictions include direct and resolved contributions, but use different 
parameters. The band in a) depicts the uncertainties which are mainly due to the long 
range matrix elements. In b) the uncertainties are indicated by the upper and lower 
curves. 
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Figure 5. Long range matrix elements as compiled in [||. The values are extracted 
from the transverse momentum distribution in pp J/^ production using different 
theoretical approximations. Only statistical errors are shown. For more details see 
text. 



The long range matrix elements are denoted as e.g. O"^ [8,^ Si] to describe the 
transition of the cc state in a colour octet configuration with angular momentum state 
^5*1. They are not calculable but are believed to be universal. Numerical values 
have been determined in a number of theoretical analyses of pp data from the CDF 
collaboration using different approximations (see |§] for a summary and references). The 
most important matrix elements for photoproduction are (9'^[8/S'o] and [8,^Pj] which 
are derived in the form of a linear combination Mk = 0^[8,^ So] + k 0^[8,^ Po]/m^, 
where is a parameter of order 3. In Figure ^, the values for Mk (upper part) and 
the matrix element (9*^ [8,'^ Si] (lower part) extracted by a few groups using different 
approximations are shown following a compilation in [^: the first group of 6 values 
is calculated in the LO collinear approximation. In the second group effects of higher 
orders have been taken into account approximately by using a parton shower Monte 
Carlo model (PYTHIA). The third group uses a kt distribution for the gluon and the 
last value in Figure ^ is derived in the kt factorisation approach. 

The theoretical band in Figure ^a is calculated for (1 < < 10) ■ 10"^ GeV^ and 
(0.3 < [8,3 Si] < 2.0) ■ 10^2 Qev^ almost covering the full range of uncertainty. The 
strong rise predicted at large z values {z > 0.7) is not seen in the data which has led to 
several theoretical attempts for an explanation. In Figure ^ calculations of a different 
group [|l^ are compared to the same data as in The authors of attempt to 



estimate effects of higher orders. Here the data are a factor 3 above the predicted sum 
of singlet and octet contributions for z < 0.7, while in Figure ^ the data are below 
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Figure 6. Left panels: Resummation of NRQCD contributions near z ~ 1 according 
to |l3| . In the upper panel predictions for pt > 1 GeV are shown and compared to 
data. In the lower panel a cut oi pt > 5 GeV is applied, right panel: Comparison of 
the cross section da/dz for pt > 2 GeV from the ZEUS Collaboration with a prediction 
from llJ using shape functions and A = 300 MeV. 



the final calculation. The rise in z due to CO contributions in Figure ^3 takes place at 
z > 0.9. This large z region is dominated by diffraction and no experimental results on 
inelastic processes are available. An intriguing question is of course a possible relation 
of such CO contributions with what is traditionally attributed to diffraction. 

The NRQCD calculations shown in Figure ^ and b neglect the energy smearing 
of the transition of the cc state to the J/tp via emission of soft gluons. This effect 
is expected to be large ai z —>■ 1. In |^^, an attempt was made to calculate this 
smearing using the technique of shape functions known from calculations in decays of 
B — > J/ip + X. The calculations for the z distributions are shown together with the 
data in Figure ^ |. In this comparison the standard cut pt > I GeV is applied. The 
smeared calculations are given for two values of A and show the expected decrease 
towards 2; — >■ 1 while the unsmeared curve (labelled "total partonic" in the figure) is 
seen to rise. However before decreasing the smeared curves are observed to increase 
above the unsmeared one (around z > 0.5). The authors of [|T^] conclude that the 

§ The authors also extract the long range transition matrix elements fom B decays. The values, which 
were used for Figure ^ agree approximately with the values derived from the pt distribution of J/ip 
mesons in pp collisions taking into account higher orders via a PYTHIA Monte Carlo simulation fTl. 
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theory is not well behaved and suggest applying a higher pt cut. Their prediction for 
Pt > 5 GeV is shown in the lower part of the figure. A comparison of this calculation 
with ZEUS results for a cut pt > 2 GeV is shown in the right panel of Figure ^ and is 
seen to give an acceptable description of the data. 



2.3. Resolved Contributions and Polarisation Analysis 

The region of low z values, z < 0.45, is explored in an analysis of HI covering a 
region 130 < W^p < 250 GeV and pt > 1 GeV. The result is shown in Figure |^a in 



comparison to predictions from . At such low z values, resolved photon contributions 
are expected and CO contributions are much larger than CS contributions. Here the 
smearing effects of soft gluon emission are not expected to play an important role. The 
number of parameters entering this comparison is larger and probably a mere cross 
section measurement will not suffice to distinguish between models. This can be judged 
from Figure |^, where the predictions within the GSM are shown with the dominant 
uncertainties, which overall amount to a factor 5. The pt distribution of the data is 
shown in Figure 0c and looks - at least with present errors - very similar to the form 
found at medium z values shown as a scaled curve. 

Measurements of J/ip polarisation in pp J/ip + X were advocated to give 
independent proof of the presence of colour octet contributions. The polarisation is 
measured via the angular distribution of the J/ ip decay leptons in the rest frame of the 
J/ip meson. Parametrising it as 1 + a cos^ a = +1(— 1) corresponds to transverse 
(longitudinal) polarisation, while a = reflects no polarisation. The data for J/ ip and 
ip' from the CDF collaboration though limited in statistics, do not show the expected 



increase of the polarisation parameter a towards higher pt predicted within NRQCDfT? 

ZEUS performed a first measurement of the J/ip polarisation in 7J9 scattering. 
Using the fiight direction of the J/ip in the laboratory system as a reference axis the 
results in Figure Ha were obtained. They are compared with the predictions of the 



"semihard" model [|T8| using unintegrated gluon density distributions based on BFKL 
evolution equations. The trend of the data is described within this model. There are 
no calculations within NRQCD/factorisation which one can compare directly with the 
data. In Figure calculated in the "recoil" system which uses the J/ip direction in 
the 7p system as a reference axis, one can however see that the a parameter increases 
with pt^^ for the Colour Singlet Model while it remains small for the octet contributions. 

Summarising, there is still no evidence of colour octet contributions in 
photoproduction at HERA. Due to the experimental and theoretical uncertainties they 
can however not be excluded. 



3. Diffractive Production of J/ip Mesons 

In the last few years it has been demonstrated that elastic J/ip production can be 
described by perturbative QCD in photoproduction and also at finite . Data for 
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Figure 7. Photoproduction of J/ij] mesons at low z from HI. Note the changed W^p 
range compared to medium z. The theoretical predictions for the z distribution in a) 
are from Kniehl et al |12| . In b) the same HI data are shown with a LO calculation 
within the CSM from |8| depicting the major uncertainties, c) da/dpf in the low z 
regime. The curve is a scaled fit to the medium z data. 




Figure 8. a) The polarisation parameter a as a function of the transverse momentum 
of the J ftp meson measured by the ZEUS collaboration . b) Theoretical predictions 
for a within the CSM (dashed curve) and the Colour Octet Model as functions of pt 
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Figure 9. Elastic photoproduction of J/ip mesons as a function of W^p (left) and 
Q"^ (right). The data from HI and ZEUS are shown with fits to the data as well as 
predictions within two pQCD models. 
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Figure 10. Elastic photoproduction of J/tp mesons: the data from HI and ZEUS are 
shown with two different theoretical predictions within pQCD using different parton 
distributions. 
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elastic diffractive processes are now available in a large kinematic range. A summary 
of the published HI data (21 pb"^) % [Tl and the new preliminary data from ZEUS 



(corresponding to ~40 pb ^) are shown in Figure ^. The data cover the regime of 
photoproduction, 0, extending almost to the kinematic limit in the photon proton 

energy W^p < 300 GeV. In the data go up to ~ 100 GeV^ 

The curves in Figure ^ show fits to the data of the form and 6 is also plotted as 
a function of in the insert of Figure The value of 6 is on average well above ~ 0.22 
expected from Regge-type fits with a soft pomeron. The distribution is fitted by a 
functional form (Q^ + M^)-", yielding n = 2.38 ± 0.11 for HI and n = 2.60 ± 0.11+°;°^ 
for the ZEUS data. 

The photoproduction data in Figure [1^ are overlaid with predictions from two 
theoretical groups using different approaches within pQCD |2l|, The basic picture 
is the exchange of two gluons between the proton and the cc pair. Apart from a number 
of technical differences the groups handle the conversion to a J/ip meson differently: 
MRT use parton-hadron duality, while FKS use a wave function. 

The important prediction is for the slope of the data which is described well by all 
calculations shown with the exception of GRV98 partons. HI has previously also found 
good agreement within the FKS calculations using CTEQ4M or MRSR2 partons [^. 

Both collaborations have extracted the effective trajectory from photoproduction 
data. The procedure is to fit the data at fixed values of t with a form The 
resulting values of a{t) are displayed in Figure |1^ with separate fits to a simple form 
a(0) + a't. The more accurate ZEUS data yield a(t) = (1.193 ± 0.01ll|^:{J}^) + (0.105 ± 
0.0241q q2o) '^^^ HI data are compatible with this result. The soft pomeron trajectory 
is ruled out. A NLO BFKL calculation of the intercept a(0)||2^ is in agreement with 
the data. 

The data shown so far are dominated by low values of \t\ since they are found 
to drop with an exponential behaviour: e~^'*' with b ^ A GeV~^. Data collected at 
high |t| are shown in Figure |l^. At high values of the proton dissociates into a low 
mass system X. ZEUS and HI, at present, measure slightly different cross sections: 
the ZEUS data are corrected to a range 0.01 < t/{t - M\) < 1 while HI corrects to 
the full kinematically allowed region of Mx- A LO calculation implementing the BFKL 
evolution equation [plf reproduces the data sets well using a value for as ~ 0.2. 



4. Summary and Outlook 

Inelastic photoproduction of J/tp mesons at HERA continues to be well described by 
the colour Singlet Model calculated in next-to-leading order. At a pt,^— 4 GeV, the 
measured cross section is reproduced by the NLO calculation, which is a factor of 
10 above the LO prediction. A reasonable alternative description of the data on a 
Monte Carlo basis can be obtained using parton distributions not integrated over the 
transverse momentum based on the CCFM evolution equations. As regards Colour Octet 
contributions, only LO calculations are available. Due to normalisation uncertainties. 
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Figure 11. The effective trajectory for photoproduction of J /ip mesons extracted 
from photoproduction data. 



ZEUS 96-97 Preliminary 




Figure 12. J/il) Production at high values of |t| in comparison with predictions from 
using BFKL dynamics. The band around the curve on the left depicts the major 
uncertainties. 
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the comparison with data at high z values is still ambiguous. NRQCD resummations at 
high z seem to necessitate a higher pt,ii) cut and then describe the data reasonably well. A 
first analysis of data at low z values shows sensitivity to resolved contributions. However 
since in this regime more parameters enter, a mere cross section measurement may not 
be sufficient to decide on Colour Octet contributions. A first attempt to measure the 
J/ip meson polarisation looks promising. This and all other analyses will profit from 
the increased statistics which will soon become available. 

For elastic J/ip production, the W^p range in photoproduction has been extended 
almost to the kinematic limit. The cross section continues to rise. Parametrising it 
as W^p, 5 is measured to be large, of order 0.8. The data at higher are in rough 
agreement with such a steep rise. Comparisons of theoretical calculations within pQCD 
show a sensitivity to the gluon distribution. A "trajectory" measurement excludes the 
soft pomeron trajectory. The measurement of diffractive processes has been extended 
to I if: I values of ~ 20 GeV^. The data is well described by calculations based on the 
BFKL equation. 
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